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ABSTRACT 


An analytical model is proposed for calculating the rotor performance and aerodynamic blade forces for 
Darrieus '^ind turbines with curved blades. The method of analysis uses a mult iple-st reamtube model, di- 
vided into two parts: one modeling the upstream half-cycle of the rotor and the other, the downstream 
half-cycle* The upwind and downwind components of the induced velocities at each level of the rotor were 
obtained using the principle of two actuator disks in tandem. 

Variation of the induced velocities in the two parts of the rotor produces larger forces in the upstream 
zone and smaller forces in the downstream zone. Comparisons of the overall rotor performance with 
previous methods and field test data show the important improvement obtained with the present model. 

The calculations were made using the computer code CARDAA developed at IREQ. The double-multiple 
st reamtube model, presented in this paper, has two major advantages: it requires a much shorter computer 
time than the three-dimensional vortex model and is more accurate than multiple-streami ube model In 
predicting the aerodynamic blade loads. 


INTRODUCTION 

The Darrieus-type vertical-axis wind turbine is 
one of the best wind energy conversion system. 
This curved- blade turbine Is now under theoreti- 
cal and experimental study in the U.S., at Sandia 
Laboratories (2-m, 5-m and 17-m turbines), for 
example, and in Canada at a larger scale with the 
230-kW machine in ..the Magdalen Islands and, 
recently a new project for building a prototype 
of about 4000-kW for the end of 1983. 

The first approach to modeling the aerodynamic 
performance of a Darrieus wind turbine was devel- 
oped by Teroplin (ref. 1), who supposed that the 
rotor is enclosed in a single s t ream tube* 

A more complex analytical model is the multiple- 
streamtube (ref. 2) In which the swept volume of 
the turbine is divided Into a series of adjacent, 
aerodynamically independent, streamtubes. The 
blade element and momentum theories are then em- 
ployed for each streamtube. This model assumes 
that the Induced velocity varies over the frontal 
disk area, both in the vertical and the horizon- 
tal directions. 

Rotor power, torque and drag are calculated by 
averaging the contributions from each stream- 
tube. The multiple streamtube model is a good 
approach for predicting the overall performance 
and axial force of Darrieus turbines under con- 
ditions where the blades are lightly loaded and 
the rotor tip speed ratios are low, (refs. 2 
and 3). 

Other types of aerodynamic performance model for 
studying vertical-axis wind turbines are based on 
the vortex theory: these are the fixed-wake 
(ref* 4) and free-vortex models (refs. 5 and 6). 
Recently, Strickland (ref. 6) extended the free- 
vortex model to the curved-blade Darrieus tur- 
bine, using a concept of the three-dimensional 
vortex structure. He used a single liftlng/llne 
vortex theory to represent an airfoil segment. 

Overall turbine performance in the form of drag, 


torque and power coefficient, have been reason- 
ably predicted by streamtube models and the ana- 
lytical results are in good agreement with exper- 
imental data* The performance values calculated 
with vortex model In the same conditions (stalled 
or unstalled blades) are not much better from 
this point of view. Although the vortex models 
have the major advantage of predicting the blade 
forces more exactly than the streamtube models 
but they have the disadvantage of requiring 
considerable computer time. 

AERODYNAMIC MODEL 

Induced Velocities 

The calculation of the induced velocities through 
the rotor is based on the principle of the two 
actuator disks in tandem at each level of the 
rotor, as shown in Figure 1. This analytical 
method uses a nuiltiple-streamtube model divided 
in two parts: one for the upstream half-cycle of 
the rotor and the other for the opposite half- 
cycle (downstream). The raultiple-streamtube 
model has been extended to the flow field 
upstream and downstream of the rotor. 

The upwind and downwind components which traverse 
each streamtube are considered separately and the 
variations in the freestream velocity are incor- 
porated into the model. The freestream velocity 
profile is given by the following relation: 

(i) 

The upwind velocity component is less than the 
local ambient wind velocity, V < and in the 
middle plane between the upstream and the down- 
stream zone there is an equilibrium-induced 
velocity, Vg<V; thus the induced velocity de- 
creases in the axial streamtube direction so that 
the downwind component is less than the equilib- 
rium velocity, V' < Vg. Figure 1 shows the ve- 
locities Induced into a pair of actuator disks in 
tandem with Lapin’s assumption, (ref. 7). For 
the upstream half-cycle of the rotor, the local 
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Lateral view Upwind view 



Figure 1 - Definition of rotor geometry for a 
Darrieus wind turbine. Two actuator 
disks in tandem. 


wind velocity decreases by an interference 

factor of u 1 and the Induced velocity becomes 

V = uV ^ (2) 

ool 

Finally the equilibrium-induced velocity is 

= (2u-l) V. (3) 

e 

For the downstream half-cycle of the rotor, is 
the Input velocity and, at the end of the stream- 
tube, the induced velocity can be written 

V* - (2u-l) V , (4) 

ool 

where u* =* VVVg is the second interference fac- 
tor for this part of the rotor, u’ < u. 

Under these conditions, the strearatube induced 
velocity is calculated by a double iteration, one 
for each part of the rotor. 

Upstream Half-Cycle of the Rotor 



Figure 2 - Angles, force and velocity vectors 
at the equatorial plane. 


asymmetrical section or a symmetrical section 
where the chord line is not tangent to the circle 
of rotation (or blade flight path), qq ^ 0, here 
we supposed oq - 0. 

Using the blade element theory and the momentum 
equation at each streamtube and by equating the 
vertical variation of the induced drag coeffi- 
cient of the rotor, we found: 

(v/v ) Rv /V ) - (V/V )] (7) 
up OQ “ L ® J 


or, In the terms of the interference factor. 


F u = Trn 
up 

(1-u) 

, n - r/R 



(8) 

where is the 

upwind conditions 

function that characterizes 

the 

F ^ 


cos9 ^ 

sln0 \ 



up 8 ttR^ 


|cos9j ‘'T 

jcos0 jcos6y 

(9) 


and 


The local relative velocity for the upstream 
half-cycle of the rotor, -7t/2 0 tt/2, is 

given by the expression: 

J(X - sin 8)^ + cos^e cos^6] (5) 

where X * mr/V represents the local tip speed 
ratio. The definitions of the angles, and the 
force and velocity vectors at the equatorial 
plane of the rotor are given in Figure 2. The 
general expression for the angle of attack is: 


a 


arcsin 


cosBcosficosa^ - (X-sin0) sin a^' 
y (X - sin6)^ + cos^O cos^6 


( 6 ) 


COSO + Cp sino 

( 10 ) 

sino - COSO 

The blade airfoil section lift and drag coeffi- 
cients, and Cq respectively, are obtained from 
NASA and Sandia Laboratories test data by inter- 
polation using both the local Reynolds number and 
the local angle of attack. 

Defining the blade Reynolds number as Re^j” Wc/o«, 
for local conditions, Re^j is given by the follow^ 
Ing relation: 

= (vc/u^) /(X - sin0)^ + cos^e cos^S (11) 


This Equation suggests the possibility of an where c is the airfoil chord and v» is the klne- 
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malic viscosity. Introducing the turbine 
Reynolds number Re^- - ioRc/Ugo, we obtain the rela- 
tionship between and Re^.: 

- ^Re^n/X^^CX^ sin6)^ + cos^0 cos^6 (12) 

For 3 given rotor geometry and rotational speed a) 
and a given velocity at each strearatube position 
from Equation (1), a value of induced local tip 
speed ratio X is chosen by assuming that the 
interference factor u is unity. Thus, Re^ and a 
will be evaluated in first approximation and the 
enarac teristlcs of the blade section profile 

Then, with Equations (10), the normal and 
tangential force coefficients are estimated, and 
Equation (9) allows us to evaluate the function 

With the first value of F^p we can calcu- 
late another value of the interference factor 
employing Equation (8) and the iterative process 
is repeated until successive sets of u are 
reasonably close. Convergence Is fast with a 
maximum number of iterations at the equatorial 
streamtube, especially for high tip speed ratios, 
and the error is less than 10”^. Once the true 
value of the Induced velocity V has been caicu* 
iated, we can obtain the local relative velocity 
W with Equation (5) and the effective angle of 
attack a with Equation (6). 

Upstream Blade Forces and Performance 

The streamwlse blade force, normal and tangential 
components may be evaluated for each streamtube 
as functions of the blade position. Half-rotor 
torque, power and drag are resolved by averaging 
the contributions from each streamtube for the 
upstream half-cycle of the rotor. 

For eacln blade in the upstream position, the non- 
dimensional force coefficients as functions of 
the azimuthal angle G are given by: 

Fj, (6) - (cH/S) (n/cosS) dc, (13) 


Thus the power coefficient for the upstream half- 
cycle of the rotor can be written: 





(18) 


Downstream Half-Cycle of the Rotor 

For the second half-part of the rotor in the 
streamflow direction (or downstream half-cycle) 
the local relative velocity is: 

^ v*^[(X* - sin6)^ F cos^ 6 cos^6] (19) 

where X’ » u)r/V 

and the angle of attack is given by Equation (6) 
where X is replaced by with m/2 $ 6 « 3m/2. 
The Induced velocity in this part of the rotor V’ 
is a function of both interference factors: up- 
stream (u) and downstream (u’), and with the con- 
dition of continuity we can obtain one of the 
velocity components in terms of_ the others for a 
certain wind velocity V , Equation (4). 


Following the same logic used for the upstream 
half-cycle, we consider that the equilibrium used 
for the velocity given by Equation (3) is the in- 
put condition of the flow in the downstream half- 
cycle at each streamtube. Thus, the iterative 
process is initialized by u' = u, where u Is the 
true value obtained for the first part of the 
rotor at each level. 


The transcendental equation which contains the 
interference downstream factor u’ becomes: 

F , u’ - 1TT1 ( 1-u’ ) (20) 

dw 

where the function F, is: 
dw 


dw 


= Nc / cose , 

8 itrJ /2 \ N |cose| T 


sln 6 X /W) 

I COS0 I cos<S / 


2 


d 0 

( 21 ) 


Downstream Blade Forces and Performance 


for the normal force, and by: 

(0) = (cH/S) (n/cosS) dc, (14) 


The normal force coefficient as a function of 0, 
for a complete blade, is : 

Fjj (0) = (cH/S) y| C’^ (w’/V„)^ (ti/cos 6) dc (22) 


for the tangential force, where C - z/H and S is 
the swept area of the rotor- 

Tlie torque produced by a blade element is calcu- 
li; ted at the center of each element. By integra- 
tion along the blade we obtain the torque on a 

complete blade as a function of 0: 

1 / 2 

T (9)=4p„cRH / C_ M'' (n/cos5) d? (15) 

up 2 “ T 

The average half-cycle of the rotor torque pro- 

duced by N/2 of the N blades is thus given by: 


and the tangential force coefficient of the blade 
as a function of 0 is given by: 

(0) = (cH/S) y* (w'/vj^ (n/cos6) d? (23) 

The torque on a complete blade, in the downstream 
half-cycle, as a function of 0 is: 

T, (0) = cRH f C' (n/cosd) dC (24) 

dw 2. 00 y , T 


^ r - 


(0) d6 


up m 

and the average torque coefficient will be: 


C 


Ql 


NcH 

2wS 


4”/' ^{hf (;fc) 


(16) 


(17) 


The average half-rotor torque has the following 
form: 



37t/2 

^dw 

nil 


(0) d8 


(25) 


and the average torque coefficient: 


1 



(26) 
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if C'X 


Q2 


The power coefficient for the downstream half- 
cycle of the rotor becomes: 


C 


P2 



^Q2 


(27) 


The power coefficient for the full cycle Is the 
weighted sum of the coefficients for the half- 
cycles: 




( 28 ) 


RESULTS AND DISCUSSION 

In order to verify this analytical model, we cal- 
culated the aerodynamic loads of the blades and 
torque as a function of the blade angular posi- 
tion as well as the power coefficient as a func- 
tion of tip speed ratio, including the effects 
due to ambient windstream shear. The theoretical 
results obtained with the double-multiple-stream- 
tube were compared with other analytical models 
such as mult iple-st reamt ube and two-dimensional 
vortex theory, and with the experimental data for 
Sandia 5-m rotor (ref. 8). The rotor of this 
turbine possesses two (or three) of the straight- 
line / circular-arc form, with a constant symmet- 
rical airfoil NACA-0015, from hub to hub. 

The double-mult iple-s Lreamtube model shows impor- 
tant retardation of the flow in the downwind zone 
of the rotor. The aerodynamic forces on the 
blades are considered as the nondimensional 
normal and tangential force coefficients for each 
element of the blade, separately, in the upwind 
and downwind positions. The elemental force 
coefficients were integrated along the length of 
the blade and the variation with the azimuthal 
angle was obtained with respect to a rotational- 
axis system, 0 yj^ where X|^ has its initial 
position linked to blade No, 1 in Figure 2. 

Figure 3 shows the resultant of the normal force 
coefficient for two blades as a function of the 
blade position. For a tip speed ratio of 
^EQ” 3*0, the zero value of Fp^ may be observed 



Figure 3 - Comparison between normal force coef- 
ficients calculated by the mulciple- 
streamtube theory, and the present 
model , 



Figure 4 - Variation with azimuthal angle 0 of 


the normal blade loading, for each 
blade, in the upwind and downwind 
zones. 



Figure 5 - Variation with azimuthal angle 9 of 

the normal force coefficient, for two 
blades, at three tip speed ratios, 

to occur at 9 <- 90° in contrast with 9 = - 90°, 
as is the case of the multlple-streamtube model. 
The variation in the normal force coefficient 
with the blade position is shown in Figure 4 for 
each blade; the normal forces are smaller in the 
downwind than in the upwind zone. The effect of 
the tip speed ratio on the normal force is pre- 
sented in Figure 5 for X^q - 1.5, 3.0 and 5.2. 





Figure 6 - Comparison between tangential force 

coefficients calculated by the inulti- 
ple-streamtube theory and the present 
model. 



Figure 7 - Variation of the tangential blade 

loading with the azimuthal angle 0 , 
for each blade, in the upwind and 
downwind zones. 



"igure 8 - Variation of the tangential force 

coefficient with the azimuthal angle 
0, for the two blades, at three tip 
speed ratios. 



Tip speed rqlio, Xfo 


Figure 9 - Power coefficient as a function of 
the equatorial tip speed ratio. 
Comparison between analytical model 
results and field test data for the 
Sandia 5-m, two-blade rotor. 



Figure 10 - Power coefficient as a function of 
the equatorial tip speed ratio. 
Comparison between analytical model 
results and field test data for the 
Sandia 5-m, three-blade rotor. 

The tangential force coefficient for two complete 
blades as a function of the azimuthal angle is 
plotted in Figure 6 for a tip speed ratio of 
XgQ ^ 3.0* The variation of this force calculat- 
ed with mult Iple-st reamtubes is similar on both 
half-cycles of the rotor. The tangential forces 
obtained with the present model are smaller on 
the downwind face of the rotor in comparison with 
the tangential forces on the upwind face. Figure 
7. The effect of the tip speed ratio on the 
tangential force is given in Figure 8 for 
XgQ = 1.5, 3.0 and 5.2. 


Figure 9 presents a comparison of the power co- 
efficient results from three analytical models: 
multlple-streamtubes , two-dimensional vortex 
theory and double-multiple-streamtubes with 
Sandia field test data for two NACA-0015 blades, 
5-m rotor. We can see that the double-multiple- 
streamtube model Is in good agreement with the 
experimental data for tip speed ratios of 3.5 to 
8.5. The results of the present model are better 
than with other streamtube methods because the 
downwind blade sees the induced velocity differ- 
ently from the upwind blade, V’<V. This con- 
cept can be improved by considering the effects 
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of the dynamic stall of the blades which Is 
prevalent at low tip speed ratios and causes the 
torque _to Increase. 

Figure 10 shows a comparison between the power 
coefficient obtained from the multiple-stream- 
Lube theory, the present model and the Sandla 
field tests for three NACA-0015 blades, 5-m ro- 
tor. There is a good agreement with test data up 
to tip speed ratio of about 7; for low speed 
ratios, X 3.5, the same phenomenon as with a 
two-blade rotor can be observed. There is a big 
difference between the mult iple-streamtube re- 
sults and experimental data, where the power 
coefficient reaches Its maximum. The tip speed 
ratio for a maximum power coefficient decreases 
mainly with Increasing drag and does not seem to 
be influenced by the number of blades. Retarda- 
tion of the flow in both upstream and downstream 
zones intensifies with Increasing rotor solidity; 
however, this effect is more significant in the 
downstream region of the rotor. 

CONCLUSION 

The present analytical model represents an exten- 
sion of the mult iple-st ream! ube model, which is 
capable of calculating the difference in the in- 
duced velocities at the upstream and downstream 
passes. This is possible If it is assumed that 
each element of the rotor is replaced by two 
actuator disks in tandem, in the wind direction. 
The aerodynamic characteristics for each element 
of the blade were therefore obtained independent- 
ly for the upwind and downwind parts of the rotor 
by using both the local Reynolds number and the 
local angle of attack. 

The double-mult iple-strearatube model allows for 
variation of the aerodynamic loads and the torque 
with the blade position for each part of the ro- 
tor. The method Is generally adequate for study- 
ing a Darrleus rotor with several curved blades. 
Field test data with a Sandia 5-m machine and 
comparisons with previous methods satisfactorily 
confirm the theoretical predictions for the over- 
all performance parameters. 

The present model could be Improved by consider- 
ing the effect of dynamic stall on the blades, 
which would enhance the ability to predict both 
loads and performance. This refinement is 
currently under way at IREQ and the results will 
be compared with those furnished by other analy- 
tical models as well as with Magdalen Islands 
field test data. 
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QUESTIONS AND ANSWERS 


I, Paraschivoiu 


Froms Art Smith 

Qi If the Darrieus produces a significant side force ^ is it not wrong to neglect the 
lateral component of velocity? 

A: No^ because the effect of the lateral component of velocity is normally smaller than 

20% of the freestream velocity (see Ref. 6); if the lateral component of velocity is 
considered the loads prediction will he better. 
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